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ABSTRACT 

■ Aims. We present new Very Large Array (VLA) radio images at 74 and 324 MHz of the SNR W44. The VLA images, obtained with 
unprecedented angular resolution and sensitivity for such low frequencies have been used in combination with existing 1442 MHz 

, radio data, Spitzer IR data, and ROSAT and Chandra X-ray data to investigate morphological and spectral continuum properties of 

^ ! this SNR. 

1 Methods. The observations were carried out with the VLA simultaneously at 74 and 324 MHz in the A and B configurations and at 

T _j_ ' 324 MHz in the C and D configurations. The radio continuum spectral index distribution was derived through direct comparison of 

| the combined data at 74, 324, and 1442 MHz. In addition, to isolate and identify different spectral components, tomographic spectral 

. analysis was performed. 

■») ■ Results. The spatially resolved spectral index study revealed that the bright filaments, both around and across the SNR, have a straight 
spectrum between 74 and 1442 MHz, with a 0.5, with two clear exceptions: a short portion of the SNR limb to the southeast, with 

^ | a varying between and +0.4 and a bright arc to the west where the spectrum breaks around 300 MHz and becomes concave down. 

We conclude that at the shell and along the internal filaments, the electrons responsible for the synchrotron emission were accelerated 

■ at the shock according to a simple diffusive shock model. The positive spectrum corresponds to a location where the SN shock is 
CIh' running into a molecular cloud and the line of sight intersects the photo dissociation region of an HII region and a young stellar object 

is present. Such spectral inversion is a classic signature of thermal absorption, either from ionized gas in the postshock region, from 
the HII region itself, or both. The curved spectrum on the westernmost bright arc is explained as the consequence of strong post-shock 
densities and enhanced magnetic fields after the interaction of the SN shock with a coincident molecular cloud. The comparison of the 
324 MHz image with a 4.5 /jm IR image obtained with Spitzer underscored an impressive correspondence between emission both to 
the north and west sides of W44, while the comparison with ROSAT and Chandra images confirm that the synchrotron radio emission 



o 



surrounds the thermal X-ray radiation. 

Key words. ISM: individual objects: W44-ISM: supernova remnants-radio continuum: ISM 

1 . Introduction cesses taking place either within SNRs (e.g. shock acceleration), 

in their immediate surrounding, or in the ISM along the line of 

The interaction of supernova remnant (SNR) shocks with inho- sight (e g therma l absorption, Lacey et al. 2001; Brogan et al. 

mogeneous interstellar gas strongly affects the spatial and tern- 2005) Based on accurate spec tral continuum index studies it is 

poral evolution of SNRs. The density distribution of the ambi- possible to constrain shock acceleration theories that predict sub- 

ent circumstellar material (CSM) and interstellar medium (ISM) tle variations in both spatially resolved (Anderson & Rudnick 

can influence the formation of a variety of SNR morphologies, m3) and inte g r ated spectra (Reynolds & Ellison 1992). In ad- 

including specific sub-classes of SNRs such as "thermal-X-ray mon> the information provi ded by low frequency observations 

composites" (also called mixed-morphology SNRs, Rho & Petre of SNRs is a powerfu i complement to higher frequency radio, 

1998), or remnants with bilateral symmetry (Gaensler 1998). x . rayi and infrared (IR) data to study the complex processes as- 

In this respect, sensitive high-resolution radio observations of sociate d with SNR shocks (Brogan et al. 2005). 
SNRs are particularly useful for studying the physics of shock 
waves traveling through nonuniform environments in the ISM. 

On the other hand, continuum radio spectra made accurate W44 (G34.7-0.4) is a bright radio SNR (S igh z ~230 Jy) with 

by extension to long wavelengths provide sensitive and often an asymmetric morphology, about (f.5 in size, and located in a 

unique means of distinguishing between different physical pro- complex region of the inner Galactic plane rich in both thermal 

and nonthermal sources. W44 is an archetype of the "mixed- 

Send offprint requests to: G. Castelletti morphology" class of remnants, characterized by a highly fil- 

* Post-Doc Fellow of CONICET, Argentina amentary radio shell and centrally concentrated thermal X-ray 

** Member of the Carrera del Investigador Cientffico of CONICET, emission. Morphological studies of this remnant in the optical 

Argentina band reveal Balmer-dominated, [Nil], and [SII] filaments from 
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radiative shocks, with some concordance with the radio emission 
(Rho et al. 1994; Giacani et al. 1997). 

This SNR is also of interest because it constitutes one of the 
few demonstrated cases of a SNR-molecular cloud interacting 
system. Seta et al. (1998) observed six giant molecular clouds 
that appear to be partially surrounding the remnant. Higher spa- 
tial resolution 12 CO observations of the full extent of W44 pre- 
sented by Seta et al. (2004) revealed that some of these clouds (at 
Vls ,r~48 km s~') are physically interacting with the remnant on 
its southeastern and western sides. These observations also show 
that the shocked molecular gas is confined to compact (~1.5 pc) 
cores located adjacent to bright radio filaments or knots. The 
physical interaction of the blast wave with the clumpy ISM is 
also supported by the detection of bright OH ( 1 720 MHz) masers 
(at LSR velocities between ~43 and ~47 km s _1 , Hoffman et al. 
2005, and references therein), and IR cooling lines from H2 
(Reach et al. 2006). 

The associated pulsar PSR B1853+01 is located inside the 
W44 shell, about 9' south of the geometric center of the rem- 
nant. A characteristic age of about 2 x 10 4 years has been esti- 
mated for this pulsar (Wolszczan et al. 1991). Frail et al. (1996) 
and Petre et al. (2002) report the detection of a radio and X-ray 
nebula powered by an associated pulsar wind. 

It has also been proposed that the EGRET source 
3 EG 1856+0114 is associated with W44 (Esposito et al. 1996). 
The 95% confidence contour for this y-ray source coincides with 
the southeastern sector of the remnant (Thompson et al. 1996). 
de Jager & Mastic hiadis (1997) explained the observed GeV y- 
rays in terms of relativistic bremsstrahlung and inverse Compton 
scattering produced by a power-law distribution of relativistic 
electrons injected by the pulsar PSR B1853+01. While this con- 
stitutes a possible explanation, a significant part of the y-ray 
flux observed may also be due to SNR shock-cloud interactions 
(Fatuzzo & Melia 2005). No y-ray radiation at TeV energies was 
detected towards W44 using either the Whipple (Leslard 1995) 
or CANGAROO (Buckley et al. 1998) telescopes, suggesting a 
spectral cutoff between GeV and TeV energies (Rowell et al. 
2000). 

A distance of about 3 kpc was first suggested for W44 by 
Radhakrishnan et al. (1972) and Caswell et al. (1975) based on 
H 1 absorption measurements. From the afore mentioned molec- 
ular studies revealing associated gas with LSR velocities be- 
tween approximately 43 and 48 km s _1 and using a flat rotation 
curve with the solar constants of Ro=7.6 kpc (Eisenhauer et al. 
2005) and Q =27.2 km s" 1 kpc" 1 (Feast & Whitelock 1997), a 
kinematical distance of 2.9 + 0.2 kpc can be derived for W44. 
Thus W44 is probably an object of the Sagittarius arm consistent 
with it being a core-collapse SNe from a massive progenitor star. 

In this paper we present the highest resolution and sensitivity 
images of the SNR W44 yet obtained at low radio frequencies. 
These images were generated from multiple-configuration Very 
Large Array (VLA) 1 observations at 74 and 324 MHz carried out 
in 2002 and 2003. Based on the comparison of the new images 
with existing VLA data at 1442 MHz, we have performed the 
first detailed study of the radio continuum spectral properties of 
W44, searching for changes in spectral index as a function of 
frequency and position. In addition, special attention was paid 
to the comparison between our high resolution radio images and 
the best existing infrared data from Spitzer and in X-rays from 
ROSAT and Chandra. 



1 The Very Large Array of the National Radio Astronomy 
Observatory is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 



2. Observations and Data reduction 

We have performed low frequency observations of the SNR W44 
at 74 and 324 MHz using multiple-configurations of the VLA. 
Table 1 summarizes various observing details. Data at both fre- 
quencies were obtained in spectral-line mode to enable radio- 
frequency interference (RFI) excision and to mitigate bandwidth 
smearing. The observations at 74 MHz and 324 MHz were per- 
formed with 64 and 32 channels, respectively. A model image of 
Cygnus A was used to correct the bandpass 2 . The absolute am- 
plitude scale were set using Cygnus A at 74 MHz and 3C 48 at 
324 MHz. Also, for the 324 MHz data, regular observations of 
the secondary calibrator 1822-096 were used for initial phase 
and amplitude calibration. At both frequencies, the data from 
each array were fully reduced and imaged separately using the 
NRAO Astronomical Image Processing Software (AIPS) pack- 
age. After initial calibration the data at 74 and 324 MHz were 
averaged down to fourteen and eight channel spectral resolution 
data cubes, respectively, before imaging. 

Two additional difficulties associated with low frequency 
radio observations are the wide-field imaging requirement im- 
posed by the large field of view (the primary beam of the 324 and 
74 MHz data are 2°5 and 1 1°5, respectively), and ionospheric- 
based phase variations that decorrelate phases on longer base- 
lines, especially at 74 MHz (Kassim et al. 1993). We there- 
fore employed wide-field imaging as implemented in the AIPS 
task IMAGR, based on a pseudo-three-dimensional multifacet 
algorithm to deal with the non-coplanarity of the visibilities 
(Cornwell & Perley 1992). Secondly, to remove ionospheric- 
induced phase distortions we employed several loops of self- 
calibration to the data from each configuration separately at 74 
and 324 MHz (Cornwell & Fomalont 1999). We note that angle- 
invariant self-calibration (as currently implimented in AIPS) is 
not sufficient to correctly compensating for ionospheric effects 
across the full field of view, especially at 74 MHz (Cotton et al. 
2004). However, self-calibration is adequate in the special case 
where the target source is located at the phase center and is the 
brightest source in the field, as is this case with W44. 

The final calibrated 324 MHz visibility data from the A, 
B, C, and D configurations were combined into a single uv 
data set and imaged using a multi-scale CLEAN algorithm with 
four different scale sizes. For the final image, all facets were 
reassembled into one large field with the task FLATN within 
AIPS, resulting in a single image with an angular resolution of 
13'.'3 x 1372, PA=-4!7, and an rms noise of 5 mJy beam" 1 . 

Data reduction and imaging at 74 MHz were performed fol- 
lowing the procedures outlined above 3 . For the A and B con- 
figuration 74 MHz data, an initial phase-only self-calibration 
was conducted using the B-configuration 324 MHz image as the 
starting model. The A and B configuration 74 MHz data sets 
were then phase-only self-calibrated; the B configuration data 
also required a final amplitude self-calibration. The individually 
calibrated A and B data sets were then concatenated, and fol- 
lowed by several more rounds of self-calibration. The final im- 
age has a synthesized beam of 397 1 x 3579, PA=-47° and an 
rms noise of 65 mJy beam -1 . 

A concern with self-calibrated images is that they can lose 
absolute position information subject to the arbitrary position of 
the starting model, which in this case may also be subject to 
ionospheric refraction. In order to check the astrometry in our 

2 See http://lwa.nrl.navy.mil/tutorial/ to obtain source models in FITS 
format. 

3 see at http://lwa.nrl.navy.mil/tutorial/ for a full description on 74 
MHz VLA data reduction. 
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Table 1. Observational summary 



Observing 


VLA 


Bandwidth 


Integration 


dates 


Config. 


(MHz) 


time (hr) 


74 MHz Parameters 


2003 Aug 31 


A 


1.562 


7.5 


2003 Sep 1 


A 


1.562 


2.6 


2002 Jun 15 


B 


1.562 


6.0 


324 MHz Parameters 


2003 Aug 31 


A 


3.125 


7.5 


2003 Sep 1 


A 


3.125 


2.6 


2002 Jun 15 


B 


3.125 


6.0 


2002 Dec 14 


C 


3.125 


4.0 


2003 Mar 16 


D 


3.125 


6.0 



74 and 324 MHz images, we compared the position of several 
small-diameter sources detected in these images with their po- 
sitions from the NRAO VLA Sky Survey (NVSS). This catalog 
has an astrometric accuracy better than 1 " in both RA and Dec 
for bright sources (Condon et al. 1998). We found and corrected 
for an average ionospheric induced refraction of 0.63 s + 0. 12 s in 
RA and 9'.' 1 + 176 in Dec at 74 MHz and 0. 14 s + 0.05 s in RA and 
3'.'52 ± 0793 in Dec at 324 MHz. 

A high resolution VLA image of W44 at 1442 MHz was pro- 
duced from interferometric data taken in 1984 and 1985 (Jones 
et al. 1993) plus short spacing information from the Bonn 100 
m 1408 MHz Survey (Reich et al. 1990). This 1442 MHz im- 
age recovers a total flux density of Si.4ghz=300 + 7 Jy and has 
a resolution and rms noise of 23" x 19" and 4 mJy beam -1 , re- 
spectively. 

To study the spectral properties of W44 (§4), special care 
was taken to match the uv coverage for the 74, 324, and 1442 
MHz VLA data. This was achieved by applying specific weight- 
ing schemes to the 324 MHz data in order to provide the 
same resolution as the 74 and 1442 MHz images, respectively. 
Similarly, we re-imaged the 1442 MHz data so that they matched 
the 74 MHz observations. Primary beam corrections were ap- 
plied to the new 74 and 324 MHz images before estimating flux 
densities. 

3. Results 

3. 1 . VLA low frequency images of W44 

Figure 1 shows our ~ 37" resolution VLA 74 MHz image of a 
portion (~ 7.5 square degrees) of the large field of view centered 
on the SNR W44. Note that while sources are detected across the 
full field-of-view, only those located near the field center are cor- 
rectly imaged by the angle-invariant self-calibration technique 
implimented here. A close-up view of the radio emission at 74 
MHz from W44 is shown in Fig. 2. This is the first sub-arcminute 
resolution image of this source ever obtained below 100 MHz. 
Though it is not sensitive to structures larger than ~ 37', it is 
sufficient to recover all information for this source about 30' in 
size. Moreover, the fact that no negative bowl is apparent around 
the SNR at this frequency confirms that most of the large-scale 
structure is recovered. However, since W44 is almost at the limit 
of the largest imaged structures, some flux might still be missing 
from the 74 MHz image. 

In Fig. 3 we present the new VLA image of W44 at 324 MHz 
with a resolution of 13" over a field of view of ~ 2.5 square 
degrees. A close-up view showing the detailed radio morphol- 
ogy of W44 at 324 MHz is shown in Fig. 4. The present image 



has a beam size 300 times better and is 80 times more sensitive 
than any previously published 330 MHz image of W44 (Kassim 
1992). The combination of all four VLA array configurations 
is sensitive to structures between 13" and 1° in size, ensuring 
complete sampling of all the spatial scales of the W44 emission. 
The high dynamic range achieved in these new observations has 
served to reveal new details previously unnoticed in W44, es- 
pecially extreme filamentary and clumpy emission that is im- 
mersed in more tenuous diffuse emission. 

At low radio frequencies the remnant has the appearance of 
a noncircular shell elongated in the southeast-northwest direc- 
tion, with an average diameter of ~30' (or 25 pc at the distance 
of 2.9 kpc adopted in this paper). The brightest emission occurs 
along the eastern boundary, a region where previous observa- 
tions have indicated the presence of molecular clouds probably 
interacting with W44. In the 324 MHz image (Fig. 4) it is appar- 
ent that the east limb is not sharp. On the contrary, faint emission 
is clearly detected at a 5<x level extending about 1' beyond the 
bright rim. A similar weak radio halo was observed along the 
eastern side of the SNR Puppis A (Castelletti et al. 2006), a di- 
rection where Puppis A is encountering a molecular cloud as is 
the case for W44. To the west, the most outstanding feature is 
the short bright arc located near 18 h 55 m 20 s , 01° 22' (J2000). 
The rest of the boundary is more diffuse in appearance. 




18 57 30 57 00 56 30 56 00 55 30 55 00 
RIGHT ASCENSION (J2000) 

Fig. 2. Radio continuum emission at 74 MHz from the SNR W44. The 
grayscale flux is linear ranging from 260 to 700 mJy bearrT 1 . The rms 
noise level is 65 mJy bearrT 1 , and the contours are traced at multiples of 
the Act noise level, 0.260 x (1, 1.76, 3.30, 4.46) Jy beanT 1 . The synthe- 
sized beam shown at the bottom left corner is 39" x 36". The white plus 
sign to the south of W44 indicates the position of the PSR B 1853+01. 

Based on the new radio observations we have estimated in- 
tegrated flux densities of 634 + 70 Jy at 74 MHz and 411 +50 
Jy at 324 MHz for W44, where the attenuation of the primary 
beam has been corrected for these estimates. The quoted errors 
take into account the statistical errors as well as uncertainties in 
the true extension of the SNR. 

In the new sensitive 324 MHz image, the low fre- 
quency counterpart of the pulsar wind nebula (PWN) around 
PSR B1853+01 discovered at higher radio frequencies and in 
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X-rays (Frail et al. 1996; Petre et al. 2002), is clearly detected 
(small white cometary feature extending to the north of the pul- 
sar position indicated by the plus sign in Fig. 4). An enlargement 
of the emission at 324 MHz associated with this PWN is shown 
in Fig. 5. Similar to the 1.4 GHz image shown by Frail et al. 
(1996), at 324 MHz the pulsar is located at the apex of an elon- 
gated structure ~ 2'5 in extent. The two maxima observed inside 
the PWN, where the intensity peaks up to 56 mJy beam -1 , were 
previously reported by Jones et al. (1993) as discrete sources. 
The integrated flux density of the PWN at 324 MHz is 0.35+0.12 
Jy. Because of the relatively lower sensitivity and angular reso- 
lution, combined with the typically flat radio spectra of all PWN, 
the W44 PWN is not apparent in the 74 MHz image. 




18 5612 11 10 09 08 07 
RIGHT ASCENSION (J2000) 



Fig. 5. A close-up image of the pulsar wind nebula around 
PSR B 1853+01 at 324 MHz. The position of the pulsar is indicated 
by the plus sign. The grayscale runs from 34 to 56 mJy beam -1 and the 
contour levels are 0.036, 0.042, 0.049, 0.056 mJy beam -1 . 



4. The radio spectral index behavior 

4.1. Integrated radio spectrum 

We have used the new estimated total flux densities at 74 and 
324 MHz together with data from the literature to determine the 
mean spectral index of W44. To construct accurate spectra of 
individual SNRs it is necessary that all flux densities be tied to 
the same absolute flux density scale. In Table 2 we list all in- 
tegrated flux density estimates for W44 between 22 and 10700 
MHz. Most of the measurements over the frequency range be- 
tween 408 and 10700 MHz were extrapolated to the Baars et 
al. (1977) absolute flux scale for consistency. In some cases an 
estimate of the correction factor was not available because the 
original reference did not list the assumed flux densities of the 
primary calibrators. Below 408 MHz (log v ~ 2.6), where the 
systematic error of the Baars et al. scale is more than ~ 5%, 
no scaling was applied. These last values have been nevertheless 
used for the spectral fitting since they do not show a large spread. 



3.5 




1 2 3 4 5 

og i/(MHz) 

Fig. 6. Integrated radio continuum spectrum of W44 obtained from the 
flux density values listed in Table 2. The filled symbols correspond to 
data from the new VLA measurements at 74 and 324 MHz, while the 
rest of the values were taken from the literature and, where possible, 
brought onto a single flux density scale. The linear fit to all of the flux 
density values yields a spectral index a=-0.37 ± 0.02. (S v oc v") and 
is shown by the line. The dashed part of the line indicates the range of 
frequencies over which no correction to bring the measured values to 
the Baars et al. (1977) scale was available. 



In Fig. 6 we show a plot of the spectrum for W44 in which 
our new flux density determinations at 74 and 324 MHz are indi- 
cated by filled circles. A single power law slope a=-0.37 + 0.02 
(S v oc v") adequately fits the flux densities measured over four 
decades in frequency without any obvious break. This integrated 
spectral index is quite flat for a typical shell type SNR, and 
agrees well with earlier estimates presented by Kassim (1989a,b, 
1992) who derived a spectral index between -0.4 < a <-0.3. 
This result confirms Kassim (1989a)'s contention that the inte- 
grated spectrum of W44 remains a power law to well below 100 
MHz. This behavior is in contrast to a number of other SNRs 
that exhibit a spectral turnover due to free-free absorption (e.g. 
W49B and 3C391, Lacey et al. 2001; Brogan et al. 2005). The 
fact that the new 74 MHz integrated flux falls nicely on the 
power law confirms that our 74 MHz image is missing little, if 
any emission on the largest spatial scales. 

4.2. Local variations in the radio spectral index 

We aim to discern whether the spectrum is a global property of 
this SNR, or if it is locally affected by the shock-ISM interaction 
and/or the presence of the pulsar in its interior. To investigate 
this, we performed a detailed study of the variations in the spec- 
trum as a function of frequency and position within the remnant. 
For this analysis we have used the new high resolution images 
at 74 and 324 MHz together with the reprocessed VLA data at 
1442 MHz. 

To accurately determine the spectral index distribution based 
on interferometric radio images, identical beam size and shape 
at the different pairs of frequencies is required. As mentioned 
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in §2, we applied appropriate t/v-tapering, in order to match the 
range of spatial scales measured at each frequency. In addition, 
to avoid any positional offsets, the images were aligned and in- 
terpolated to identical projections (field center, pixel separation, 
etc). The final common resolution is 40" x 36", PA=-62° for 
the 74/324 MHz comparison and 14" x 14", PA=-45° for the 
324/1442 MHz pair. 

The spectral study was carried out by direct comparison of 
the matched images at the three frequencies and by construct- 
ing tomographic images. The different procedures ensure that 
the observed spectral features are robust, and unaffected by any 
artifacts or biases in the databases used for the comparison, like 
zero level differences, smooth background variation, etc. 

In spite of the relatively good fidelity attained in the new im- 
ages, some residual small-scale fluctuations, probably not real, 
were apparent at the highest spatial resolution traced in the spec- 
tral index maps. We therefore preferred to further smooth all 
three frequency maps to a 50" resolution, thus avoiding any 
small scale variation that could mask the main spectral features. 
We based our analysis on spectral comparisons derived from 
these smoothed databases. 

4.2.1. Spectral index maps 

Figure 7a shows a 324 MHz continuum image of the SNR W44 
smoothed to 25" resolution for comparison, while Fig. 7b and 
Fig. 7c display the spatial spectral index distribution in the rem- 
nant obtained from the direct ratio of the matched 324/1442 
MHz and 74/324 MHz data, at the angular resolution of 50". 
To construct these spectral index maps the images were clipped 
at the 4cr level of their respective noise levels. Since, as noted 
in §4.1, the global spectrum of W44 remains a power law to 
well below 100 MHz, we use the same scale in Fig. 7 for both 
frequency pairs in order to easily identify continuum spectral 
changes that depend on both position and frequency (spectral 
curvature). 

In Fig. 7b it is readily apparent that the brightest radio 
emission all along the periphery of W44 has a spectral index 
°324 2 ~ ~0.4, -0.5, compatible with a first-order Fermi mech- 
anism operating at the shock. In addition, the bright filaments 
observed across the SNR have a similar a, in contrast with the 
slightly flatter diffuse interior. 

In the spectral index map traced between 74 and 324 MHz 
(Fig. 7c) the most striking feature is the short band seen on the 
southeastern limb (near 18 h 56 m 50 s , 01° 17'), where the spec- 
tral index becomes positive, with a varying between and +0.4. 
This spectral inversion occurs at a site where previous studies 
have demonstrated that the SNR is physically interacting with 
a dense molecular cloud (Seta et al. 2004; Reach et al. 2005). 
Also, as shown below (§5.1) from 24 fim and 8 fim Spitzer ob- 
servations, this is the location of an HII region (named G034.7- 
00.6 in the HII regions catalog by Paladini et al. 2003, cen- 
tered at 18 h 56 m 47.9 s , 01° 17' 54" J2000) with the IRAS point 
source 18544+01 12 (identified as a young stellar object) located 
on its border. Thus, the spectral inversion can be interpreted 
as originating from low frequency free-free absorption at the 
SNR/molecular cloud shock boundary. Absorbing ionized gas 
could have been created by the passage of a fast, ionizing J-type 
shock that passed through the molecular cloud and dissociated 
and ionized the gas. A similar spectral inversion was observed 
by Brogan et al. (2005) in the SNR 3C 391, where a spectral 
turnover for frequencies below 100 MHz was found at the posi- 
tions where the SN shock was running into dense molecular gas. 
An alternative explanation is that thermal absorption is occuring 



inside the boundaries of the coincident HII region, along its pe- 
riphery where the thermal electron density might be highest. We 
note that this localized absorption has a negligible influence on 
the total integrated flux and thus has no measurable impact on 
the integrated continuum spectrum. 

From the comparison between Fig. 7b (324/1442 MHz) and 
Fig. 7c (74/324 MHz) one could have the misleading impression 
that the spectral behavior is reversed from one frequency range 
to the other. In fact the spectral index of the limb and bright 
filaments is the same (near ~ -0.4 to -0.5) in both frequency 
pairs, and it is the diffuse interior that changes to a steeper value 
in the 74/324 MHz comparison. This effect, however, needs to 
be confirmed on the basis of more sensitive 74 MHz data before 
advancing an interpretation. 

What is readily apparent from the comparison of Fig. 7b and 
Fig. 7c is that the westernmost arc changes its spectrum from 
a yiA — 0-4 to a flatter value a^ 4 ~ -0.2. This feature is also lo- 
cated at a site where the SN shock is probably running into a 
dense molecular cloud, where the intense near infrared emission 
(see §5.1) and the bright optical filaments observed in the [SII] 
and [Nil] lines (Mavromatakis et al. 2003) suggest the presence 
of radiative shocks. In this case, it is possible that the observed 
change in the radio continuum spectrum is due to absorption by a 
layer of thermal electrons associated with the sheets responsible 
for the filaments. 

4.2.2. Tomographic images 

To investigate the fidelity of our radio continuum spectral index 
maps, we also produced a tomographic gallery (see for example 
Katz-Stone & Rudnick 1997). A tomographic map between two 
images at different frequencies is made by using a test spectral 
index a t to scale the brightness in the higher frequency image 
(5 Vl ), and then subtracting this scaled image from the lower fre- 
quency image (S v ,). Thus, a tomographic gallery is obtained by 
calculating S t = S Vl - (v2lv\) ai S Vl , stepping through a range 
of test spectral indices a t . Features which have a spectral index 
identical to the test value will disappear in the S, image. Spatial 
components that have different spectral indices will appear as 
positive or negative residuals depending upon whether the spec- 
trum is steeper or flatter than the assumed a t value. Such a pro- 
cedure allows us to trace fine-scale spectral index changes, and 
is particularly useful for disentangling the superposition of com- 
ponents with different spectra that can overlap along the line of 
sight. 

In Fig. 8 we present a tomographic gallery of 324 vs. 1442 
MHz data, while the comparison between 74 and 324 MHz data 
is shown in Fig. 9. In both cases the tomographic images were 
generated for five test spectral index values: a t =+0.1, or t =-0.1, 
ff t =-0.3, a' t =-0.5, and a t =-0.1. In these maps the positive spa- 
tial components (seen as light regions) are steeper in spectral in- 
dex than the assumed test value, while the negative (dark) residu- 
als correspond to features having a spectrum flatter than the test 
value. In both figures, the bottom right frame displays the 324 
MHz continuum image for comparison. 

In the 74/324 MHz comparison (Fig. 9), the feature with pos- 
itive spectrum noted earlier near ~ 18 h 57 m , 01° 17' is strik- 
ingly revealed as a short, dark filament in the eastern border 
of all the tomographic images. It is especially evident in pan- 
els from fft=-0.1 to ^=-0.3 as a structure with flatter spec- 
trum than these values. Also, the tomographies traced between 
74 and 324 MHz emphasize the westernmost arc, that stands out 
as a flatter spectral index feature in the images traced at a t =-0.5 
and a t =-0.7 and disappears against the background at a t =-0.3, 
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while for the tomographies traced between 324 and 1442 MHz 
(Fig. 8) the best fit is attained for a t between -0.5 and -0.7, 
thus confirming the negative (concave down) curvature of the 
spectrum recognized earlier. 

In summary, all of the important radio continuum spectral 
index characteristics discussed in the previous section are con- 
firmed through this alternative analysis. 

An important conclusion from this new radio continuum 
spectral study is that neither spectral index maps nor tomogra- 
phies reveal any trace of coupling between the pulsar and the 
surrounding plasma beyond the PWN; that is, there is no spectral 
signature, such as a smooth radial steepening, that would suggest 
a progressive aging of particles from their injection near the pul- 
sar to their transport further out into the SNR shell. Moreover, 
the spectral index obtained for the shell (a — 0.5), compatible 
with the diffusive shock acceleration mechanism in the test par- 
ticle limit, confirms that the relativistic electrons responsible for 
the shell radio emission are accelerated right at the shock front, 
with no influence from the central neutron star. Also, no spec- 
tral signature is observed that can be related to the presence of 
the 95% confidence ellipse corresponding to the EGRET y-ray 
source that covers most of the southeast corner of W44. 

5. Multi-wavelength Comparisons 

5. 1 . Radio, infrared and molecular emission 

Recently, Reach et al. (2006) used Spitzer Space Telescope 
GLIMPSE data to resolve the mid-infrared structure of W44 be- 
tween 3.6 and 8 pm (Benjamin et al. 2003). To study the rela- 
tionship between the infrared and radio features in the SNR W44 
we have combined the data at 4.5 pm with the new radio image 
at 324 MHz in a false color image (Fig. 10). In this figure, the 
mid-infrared emission is shown in orange and the radio emis- 
sion in blue; magenta thus corresponds to regions where both 
spectral bands overlap. The distribution of compact OH (1720 
MHz) maser spots (Claussen et al. 1997; Hoffman et al. 2005) 
are shown as plus signs. 

The mid-IR (MIR) image reveals extensive filamentary emis- 
sion and some isolated bright clumps. This emission is particu- 
larly bright in the north and southwestern parts of W44, where 
bright near infrared H2 line (2.12 pm) emission is also present 
(Reach et al. 2005). Figure 10 reveals the impressive agreement 
between the low frequency radio emission and the infrared 4.5 
pm emission to the north and southwest, where the brightest 
radio filaments are almost exactly coincident with the infrared 
radiation. In contrast, the central and eastern regions do not 
show much obvious correspondence with the exception of a few 
bright knots. The westernmost bright radio arc near 18 h 55 m 20 s 
01° 22', which is also bright in U a , [SII] and [Nil] lines (Giacani 
et al. 1997; Mavromatakis et al. 2003), has an infrared counter- 
part that mimics the appearance of the radio arc. 

The mid-IR study is also useful for understanding W44 in the 
context of the surrounding interstellar medium. Figure 1 1 dis- 
plays a three color image of W44 comparing the radio morphol- 
ogy at 324 MHz (blue) with Spitzer GLIMPSE and MIPSGAL 
images at 8 pm (green) and 24 pm (red), respectively (Carey et 
al., in preparation). These infrared bands are excellent tracers of 
the warm dust associated with star formation. The most conspic- 
uous feature in Fig. 1 1 is the ~ 10 arcmin region just outside the 
southeast limb of W44 (red nebula), identified by Paladini et al. 
(2003) as the HII region G034.7-00.6. In this case the combina- 
tion of the 8 and 24 pm images reveals with impressive detail the 
hot dust grains (in red) in the Stromgren sphere, limited to the 




Fig. 10. A color composite image showing the spatial correlation be- 
tween the mid-infrared emission as observed by the Spitzer Space 
Telescope at 4.5 pm (in orange) from Reach et al. (2006), and the new 
low frequency radio image at 324 MHz (in blue). Features where both 
spectral bands overlap are magenta in color. The positions where OH 
(1720 MHz) maser emission was detected are indicated with white + 
symbols (Hoffman et al. 2005). 



east by an annular photo dissociation region (PDR) dominated 
by polycyclic aromatic hydrocarbon (PAH) emitting near 8 pm 
(green). The bright yellow spot abuts the only portion of the east- 
ern limb of W44 with a concave shape (near ~ 18 h 57 m , 01° 17') 
coincident with the brightest radio filament and the region de- 
termined to have a positive (inverted) radio spectrum (as shown 
in §4.2). The IR point source is the IRAS source 18544+0112, 
proposed by Molinari et al. (1996) to be a young stellar object 
(YSO) in a primitive evolutionary stage. 

As summarized in the introduction (§1), W44 is one of the 
few SNRs known to be physically interacting with its parent 
molecular cloud complex, with high-angular resolution molec- 
ular studies in different near-IR (NIR) and millimetric lines (see 
for example Seta et al. 2004; Reach et al. 2005) clearly indicat- 
ing its evolution within a giant molecular cloud. In particular, 
one molecular cloud component has been detected adjacent and 
closely parallel to its bright eastern radio boundary (Reach et al. 
2005, see for instance Fig. 3e in their paper), while another one is 
located near the very bright western arc. Seta et al. (2004) pro- 
vide a useful schematic illustration of the various CO features 
discovered towards W44 (Fig. 3 in their paper). In particular, 
and near the radial velocity of ~ 47 km s _1 approximately corre- 
sponding to the 2.9 kpc distance attributed to W44, they identify 
a CO "edge" feature exactly at the position of the YSO and the 
indentation of the eastern rim (E2 in their nomenclature). Based 
on the abrupt changes observed in the physical parameters, Seta 
et al. (2004) interpret the SNR as interacting with this particular 
molecular cloud component. 

In this scenario involving molecular gas, an SNR, an HII re- 
gion, and a YSO, it is worthwile to ask whether this may be a 
case of star formation triggered by the action of a SN shock on 
a coincident molecular cloud or, alternatively, whether it might 
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illustrate a case of sequential star formation at the periphery of 
the HII region. This last case has been, for example, observed 
in the Sh 217 and Sh 219 HII regions (Deharveng et al. 2003). 
Although Molinari et al. (1996) proposed a kinematic distance 
of 3.7 kpc for the YSO, this distance might be an overestimate 
and the YSO is instead co-distant with W44, at about 2.9 kpc. 
The question then is whether the age of W44, of a few 10 4 yrs, is 
sufficient to have been responsible for the formation of new, mas- 
sive stars. The timescale of star formation is currently under am- 
ple debate, but it is generally agreed that a few 10 6 yrs are nec- 
essary for the formation of protostars from condensed molecular 
matter (e.g. Tassis & Mouschovias 2004, and references therein). 
We are therefore led to the conclusion that if the YSO is related 
to W44, then the cloud collapse must have been initiated by the 
action of the precursor star's wind and not directly by the SNR. 
Therefore we feel it is more plausible that the YSO is associ- 
ated with the PDR around the HII region G034.7-00.6. Detailed 
NIR, MIR and millimeter studies are planned to further investi- 
gate the star formation in this complex region, with the goal of 
uncovering the underlying genetic links between the SN shock, 
the HII region, and the protostar. 

From Fig. 11 regions of more nebulous infrared emission 
toward the north-northwestern sides of W44 are also apparent, 
and the yellow spot about ~ 1 arcmin in size seen to the east of its 
geometrical center corresponds to the HII region G034.7-00.5 
(centered at 18 h 56 m 26.5 s , 01° 20' 38", Paladini et al. 2003). 
Observations of radio recombination lines (RRL) would clearly 
be helpful to further elucidate the presence of thermally emitting 
gas in this area. 

5.2. Radio and X-ray emission 

The X-ray emission from W44 has been investigated using 
the Einstein, EXOSAT, ROSAT, ASCA, and Chandra telescopes 
(Smith et al. 1985; Jones et al. 1993; Rho et al. 1994; Harms 
et al. 1996; Shelton et al. 2004). In Fig. 12a we present a com- 
parison between the 324 MHz radio image and the X-ray emis- 
sion as observed with ROSAT in the energy range 0.2 to 2.4 
keV (taken from ROSAT archives). As reported by Rho & Petre 
(1998), the brightest X-ray features are preferentially located in 
the center of the remnant, in a region of relatively low radio sur- 
face brightness. As expected, no correspondence is observed be- 
tween the mostly thermal X-ray emission (Shelton et al. 2004) 
and the filamentary synchrotron radio emission. A number of 
physical processess such as cloud evaporation, electron thermal 
conduction, and entropy mixing have been proposed to explain 
the observed X-ray properties (see Shelton et al. 2004, for a re- 
view of the models). Reach et al. (2005) conclude that probably 
both, evaporating clumps and thermal conduction behind radia- 
tive shocks, are operating to produce the observed X-ray radia- 
tion. 

Fig. 12b shows the Chandra image of W44 produced by 
Shelton et al. (2004) from the combination of 6 different point- 
ings with the outer contour (traced at Act level) of the 324 MHz 
radio emission plotted for comparison. Shelton et al. (2004) sug- 
gest that diffuse X-ray emission is observed extending slightly 
beyond the northwestern border of the radio synchrotron emis- 
sion. The new high-sensitivity radio images presented here, how- 
ever, do not support this conclusion. 

6. Energetics of SNR W44 

As W44 harbours a pulsar with a PWN and is spatially coin- 
cident with the GeV gamma-ray source 3EG 1856+0114, it is 



useful to calculate its total energy content and compare it with 
other Galactic SNRs with central compact objects. The total syn- 
chrotron energy is stored in the magnetic field as well as in the 
kinetic energy of the relativistic particles. Under the assumption 
of equipartition between particles and magnetic energy, the min- 
imum energy content and the corresponding magnetic field can 
be obtained from the following relations (Moffet 1975) 

U min =0.50(flAL) 4/7 V- 3/7 (1) 

B min = 23(aAL/V) 2/1 (2) 

where V is the volume of the source, A is a factor which 
depends on the assumed lower and upper cutoff frequencies and 
the spectral index, and L is the total luminosity of the source. We 
note that electrons may not be the only energetic particles within 
the source; an appreciable amount of energy may also be stored 
in energetic baryons. The total particle energy can be described 
as Up— a U e , where a, the energy ratio between electrons and 
baryons is usually assumed to be ~50 (Moffet 1975). 

Using the integrated spectral index of -0.37 and assuming 
isotropic radiation, we derive a total luminosity for W44 of 7.8 
x 10 34 erg s . The volume of the source is estimated to be 
about 2 x 10 59 cm 3 by approximating W44 as a sphere of ra- 
dius 11.6 pc. For typical lower and upper cutoff frequencies of 
10 7 and 10 11 MHz, the total minimum energy content is 5.8 x 
10 49 erg. Following the same procedure and with similar as- 
sumptions we also re-estimated the total energy content of the 
PWN based on the 324 MHz data. We derive its minimum en- 
ergy as U,„i„(PWN) = 1.2 x 10 47 erg, in good agreement with 
Frail et al. (1996)'s results (£/„„■„ - 1 x 10 47 erg) estimated from 
higher radio frequencies observations. 

The minimum energy derived for W44 is comparable to the 
value obtained for SNRJCTB80 (U mi „ = 5 x 10 49 erg, Castelletti 
et al. 2003), another SNR containing a pulsar, and two orders 
of magnitude higher than the energy content of W28 (t/ m ,„ = 
2 x 10 47 erg, Dubner et al. 2000). 

The corresponding minimum magnetic field in the SNR is 
B ml „=l 8 //G, consistent with the compression of a typical ISM 
magnetic field (about 5 yuG, Wielebinski & Beck 2005) by a fac- 
tor of three to four. 

7. Conclusions 

We have presented new high-sensitivity and high-resolution low 
frequency VLA radio images of the SNR W44. These images 
provide an unprecedented view of W44 with an angular reso- 
lution of ~37" and ~13" at 74 and 324 MHz, respectively, an 
improvement of over two orders of magnitude over the best, pre- 
vious low-frequency image of this SNR. 

Based on our 74 MHz image sensitive to structures with an- 
gular scales from a few arcseconds to ~ half a degree, we mea- 
sure an integrated flux density for W44 of S74MHZ -634 Jy. 
Similarly, from our 324 MHz image sensitive to spatial struc- 
tures between ~13 arcsec and ~1 degree, we derive an integrated 
flux density S 324 mhz -411 Jy. We have used these measure- 
ments, together with carefully considered flux densities present 
in the literature, to calculate an accurate global radio spectrum 
of the SNR. The best fit to the data yields an integrated spectral 
index ff=-0.37+0.02, constant over four decades in frequency 
with no significant turnover down to 74 MHz. The straight power 
law spectrum is in contrast to a number of other SNRs that ex- 
hibit spectral turnovers at low frequencies, and is consistent with 
previous low frequency measurements of this object. 
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The new 74 and 324 MHz data, together with existing 1442 
MHz data, were combined to perform the first careful, spa- 
tially resolved study of the synchrotron spectrum in W44. High- 
dynamic range observations at three well spaced frequencies are 
ideal for accurately constraining the variation of the spectrum 
within the source. Our results were obtained first from standard 
analysis of spectral index maps, and thereafter confirmed using 
spectral tomography techniques. The main findings of our radio 
continuum spectral index study are: 

(i) The eastern limb, which is the brightest portion of W44 
SNR, as well as the narrow filaments observed in total power 
across most of the surface of the source have a straight power 
law spectrum between 74 and 1442 MHz with a — 0.4, -0.5. 
This result is consistent with the predictions of a simple, diffu- 
sive shock acceleration model applied to the filamentary shell of 
W44. 

(ii) From the 74/324 MHz comparison we have discovered a 
narrow feature with a positive spectrum (0.0 < a < +0.4 ) located 
towards the southeast border of W44, coincident with an inden- 
tation in the limb of the SNR. At this location existing studies 
have shown that the SNR shock is interacting with an external 
molecular cloud. Also, flanking this portion of W44, mid-IR ob- 
servations have revealed the presence of an almost circular HII 
region surrounded by a PDR region. Exactly at the interface be- 
tween this HII region and W44, as seen in the plane of the sky, 
a young stellar object is present. A likely explanation of the in- 
version in the low frequency radio continuum spectrum in this 
region is free-free absorption from ionized gas in the post shock 
region at the SNR/molecular cloud interface. If this is the case, it 
marks the second case of spatially resolved low frequency radio 
observations tracing an SNR/molecular cloud interaction, fol- 
lowing Brogan et al. (2005)'s detection of thermal absorption 
at the periphery of 3C 391. An alternative explanation is that the 
thermal absorption is occuring at the borders of the cavity carved 
by the collindant HII region, where the column density of ther- 
mal electrons might be highest. The spectral inversion is prob- 
ably produced by a combination of both these effects. However 
this localized absorption is negligible compared to the total in- 
tegrated flux and thus has no measurable influence on the inte- 
grated continuum spectrum. Further IR and millimeter observa- 
tions are planned to investigate the 3-D distribution of the gas 
emitting in the different spectral regimes, and to elucidate the 
precise physical processes acting at this site. This should shed 
light on the complex chain of interconnected events that can tie 
together the evolution of W44, its parent molecular cloud, pre- 
cursor star, and neighbouring HII region in a scenario consistent 
with the formation of a new generation of stars exactly at the 
interface SNR/molecular cloud. 

(iii) The comparison of the spectrum between 74 and 324 
MHz and between 324 and 1442 MHz revealed a spectral flat- 
tening ocurring on the westernmost arc of W44 at the lowest 
frequency pair (concave down spectrum). Our hypothesis is that 
this spectral flattening is related to the SNR shock colliding with 
a molecular cloud in this region, consistent with the bright op- 
tical filaments and IR emission observed there. In this scenario 
the spectral flattening could be a consequence of strong post- 
shock densities and enhanced magnetic fields. If the electron en- 
ergy spectrum hardens at lower frequencies, regions in which the 
magnetic field is higher than the surroundings will appear both 
brighter and with a flatter spectrum. A firm estimate of the dis- 
tribution of magnetic field strength throughout the SNR would 
be important to test this hypothesis. 

(iv) The present study revealed that there is no evidence in 
the radio continuum spectrum of any coupling between the as- 



sociated pulsar PSR B1853+01 and the surrounding SNR shell 
(that could, for example, have been observable as a gradual 
steepening from the pulsar to the shell). Furthermore, there is 
no spectral signature indicating any link to the EGRET y-ray 
source whose 95% confidence ellipse covers most of the south- 
east corner of W44. 

We have also used the new 324 MHz image to perform 
the first detailed low frequency radio/IR comparison for W44. 
Bright filamentary IR emission detected by Spitzer at 4.5 //m 
shows an impressive correspondence with the radio emission to- 
wards the north and western portions of the remnant. To the cen- 
tral and eastern regions of W44, however, almost no IR counter- 
part is observed for the radio emission. 

From the comparison of the new radio image at 324 MHz 
with the X-ray emission as seen by ROSAT between 0.2 and 2.4 
keV and with a high-resolution Chandra image, we can con- 
firm that the synchrotron radio emission appears to surround the 
thermal X-ray plasma. As determined from previous studies, the 
X-ray emission is probably related to a combination of the evap- 
oration of clouds and thermal conduction occuring behind ra- 
diative shocks. In addition heavy absorption likely produced by 
foreground molecular clouds may explain the absence of X-ray 
radiation towards the southeast. 

We have also investigated the energetics in W44 based on the 
new observations and accurate, re-derived integrated spectrum. 
We find that the total minimum energy content is 5.8xl0 49 erg, 
a value in agreement with the estimates for other SNRs contain- 
ing central pulsars. The minimum, equipartition magnetic field 
estimate is 18 yuG. 

Finally, we note that this study adds to a growing list of ex- 
amples confirming the importance of high resolution, low fre- 
quency radio observations in providing an important new tool 
to demonstrate and understand the physical interaction between 
SNR shocks and the surrounding dense material. This also points 
to the great potential of an emerging generation of much more 
sensitive low frequency instruments. 

Future work will involve a detailed comparison of the radio 
spectral index distribution and the molecular emission in order to 
refine the analysis of the interaction between W44 and its parent 
molecular cloud. New IR and optical studies are also planned 
in order to improve our understanding of the underlying genetic 
links and causal relationships between the evolution of the SNR 
and the adjacent YSO in this inner Galactic complex. 
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RIGHT ASCENSION (J2000) 



Fig. 1. This image shows a portion of the primary beam of the VLA at 74 MHz centered on the SNR W44, after the combination of A- and B- 
configuration data. The image displayed does not include primary beam correction. The field covers ~7.5 square degrees at an angular resolution 
of 39" x 36" (P.A.=-47°). The rms noise is 65 mjy beam -1 . 
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RIGHT ASCENSION (J2000) 

Fig. 3. An image of the field of view around the SNR W44 at 324 MHz. The image displayed does not include primary beam correction. The 
grayscale varies linearly between 15 and 70 mjy beam -1 . The angular resolution is about 13" x 13" with a position angle of -4°7 and the 
sensitivity level is 5.0 mjy beam -1 . The source located near the bottom right corner (~ 18 52 30, +00 20) is the SNR Kes 79. 
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RIGHT ASCENSION (J2000) 



Fig. 4. Image of SNR W44 at 324 MHz. The grayscale is linear ranging from 20 to 60 mjy beam , The rms noise level is 5.0 mjy beam '. The 
synthesized beam shown at the bottom left corner is 13" x 13". The black plus sign indicates the position of the PSR B 1853+01. 




Fig. 7. (a) 324 MHz image of W44 with 25" resolution showing the location of the brightest filaments, (b) Spectral index map between 324 and 
1442 MHz (50" resolution), (c) Spectral index map between 74 and 324 MHz (50" resolution). The 0. 15 Jy beam -1 contour from the 25" resolution 
324 MHz image is included on each panel to facilitate the comparison between spectral continuum and total power features. Only regions with flux 
densities greater than 4<x were used to create the spectral index maps in (b) and (c). Both spectral index maps have the same color scale (displayed 
to the right). 
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Table 2. Integrated flux densities on the SNR W44 
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No correction to Baars et al, (1977) scale was applied. 

Flux density scale from VLA Calibrator Manual, http:/www.aoc. nrao.edu/~gtaylor/calib.html. 
The correction factor was not available. 
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Fig. 8. Gallery of tomographic images for the SNR W44 constructed from data at 324 and 1442 MHz, which were matched in uv coverage. The 
data at both frequencies have a common 50" resolution. The spectral index a t is indicated at the top left corner of the panels. Bright areas have 
radio continuum spectral indices steeper than the assumed a t and dark areas imply flatter spectrum. The grayscale is the same in all panels. A 
324 MHz intensity contour at 260 mJy beam -1 is plotted for reference. The image of W44 at 324 MHz plotted for comparison is displayed in the 
bottom right frame. 
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Fig. 9. The same as Fig. 8 but between 74 and 324 MHz. Bright areas have radio spectral indices steeper than the assumed a t and dark areas imply 
flatter spectral indices. The greyscale is the same for each spectral index panel. The 324 MHz contour at 140 mJy beam -1 is also plotted in each 
panel for reference. A greyscale 324 MHz image with a range of 26 to 70 mJy beam -1 is shown in the lower right panel for comparison. 




Fig. 11. A high resolution comparison of the radio continuum emission at 324 MHz of the W44 and the infrared emission in the SNR region. The 
green and red images corresponds to the Spitzer 8 and 24 jim data, while in blue the 324 MHz image of the remnant with 13" resolution is shown. 
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Fig. 12. An X-ray/radio comparison of the SNR W44. (a) The grayscale corresponds to the thermal X-ray emission observed with ROSAT in the 
0.2-2.4 KeV range, while the overlaid contours trace 324 MHz radio emission at 20 and 50 mjy beam -1 , (b) Chandra 0.7-2.6 keV image of W44 
(as taken from Shelton et al. 2004) with the 20 mjy beam -1 (4<x) radio CONTOUR at 324 MHz superimposed for comparison. The white plus sign 
shows the position of the pulsar PSR B 1853+01. 



